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Abstract — In this paper, we present a novel and 
computationally efficient water filling algorithm for multiuser in 
cognitive radio network (CRN). This algorithm is based on the 
multiuser water filling theorem and persuades the subcarrier 
allocation for a multiple access system. This approach 
exaggerates the total bit rate under the constraints of 
user-individual power budgets. Acquiring Energy Efficiency 
(EE) is very troublesome for wireless communication systems 
exclusively for CRN. The analysis here is with the mean EE 
maximization problem for secondary user (SU) with both 
primary user (PU) and SUs. The solution to the problem starts 
by examining both the peak and mean transmit power 
constraints for the SUs and outage probability constraint for the 
PUs. The problem here is nontrivial. The results attest that the 
efficiency of the system is enhanced with the proposed water 
filling algorithm and also it has been observed that the outage 
probability is reduced. Further, there is enhancement in the 
capacity of a MIMO system. 

Index Terms — Cognitive radio networks, energy efficiency, 
outage constraint, power budget, water filling algorithm 


I. INTRODUCTION 

Cognitive Radio (CR) is an adaptive, intelligent radio and 
network technology that can automatically detect available 
channels in a wireless spectrum and change transmission 
parameters enabling more communications to run 
concurrently and also improve radio operating behavior. It is 
a form of wireless communication in which a transceiver can 
intelligently detect which communication channels are in use 
and which are not, and instantly move into vacant channels 
while avoiding occupied ones. This optimizes the use of 
available radio-frequency (RF) spectrum while minimizing 
interference to other users. This technology obviously 
promotes the spectrum efficiency (SE) among primary user 
(PU) and secondary user (SU) by spectrum sharing. With the 
condition that SUs interference to PUs is restricted so that the 
licensed user’s performance depravation is in tolerable 
range[l] and so the transmission between PUs and SUs occurs 
simultaneously. The EE grew as a matter of point for future 
wireless networks since the cost of energy is promulgated on 
both financial and ecological burden on its development[2] 
beyond the SE unfortunately. Increase in EE is a fundamental 
necessity for SUs in order to accomplish improved utilization 
of limited power transmission which consumes to maximum 
SE and spectrum sensing. 
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Some more works are designed for a near optimal algorithm 
with low complexity to maximize the SUs EE under some 
practical constraints [3], [4]. As it is known that 
Multiple-Input Multiple-Output (MIMO) systems are used to 
get higher data rate as compared to a normal SISO system, the 
same power budget and signal to noise (SNR) is considered. 
A comparison of MIMO system with a SIMO reveals that the 
MIMO system need lesser transmit power than the SISO 
system in order to achieve the same capacity. As we need to 
minimize the energy consumed by the circuit and want to 
maximize the capacity of a system and that is possible only if 
we use multiple MIMO system. 

Water filling algorithm is proposed to maximize the SUs EE. 
The EE design for slow fading scenarios is suitable to perform 
static optimization depending upon the instantaneous channel 
state information (CSI). It is impossible to provide some 
performance for each state in fast fading channels to employ 
statistical metrics such as SU’s mean EE and PUs outage 
probability constraint. The capacity of a MIMO system can 
further be increased if we know the channel parameters both 
at the transmitter and at the receiver and assign extra power at 
the transmitter by allocating the power according to the water 
filling algorithms to all the channels. 

The concentration is on SE with SU’s capacity in fast fading 
channels. Maximizing SU’s capacity and outage probability 
of PU’s with statistical CSI do not help in increasing EE 
directly[5]. Here, we maximize SU’s mean EE in fast fading 
channels where SUs and PUs are involved and also we 
propose a water filling algorithm to address the optimal power 
allocation problem and to achieve the parameters like 
bandwidth, delay, jitter and throughput of Quality of Service 
(QoS)[6]. To safeguard the QoS of PU’s transmission, we 
consider PUs outage probability constraint. Also we consider 
peak and mean transmit power constraints of SUs. The 
problem identified here is formulated in fractional 
programming initially. It is observed to be nontrivial since the 
objective function is non convex[7]. 

To solve it, we initiate transformation of problem to a 
parameterized convex optimization problem based on 
fractional programming. And by using Lagrangian duality 
theory, the problem is decomposed into multiple sub 
problems in parallel. Then an efficient water filling algorithm 
has been proposed to derive optimal power allocation strategy 
for SUs in order to maximize its mean EE. This maximizes the 
total bit rate under the constraint of a maximum transmit 
power per user. An algorithm which computes the power 
spectral density for each user has been presented recently 

[14]. 

II. SYSTEM MODEL 

In underlay CRN, the secondary transmission link coexists 
with primary transmission link is considered. Here, PT stands 
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for PU transmitter, PR for PU receiver, ST and SR for SU 
transmitter and SU receiver respectively. Because of its 
spatial reuse property, we consider underlay CRNs. The links 
are ergodic, stationary and mutually independent in fast 
fading channels between the transmitters and receivers of PUs 
and SUs. 



ST 


SR 


Fig. 1 System model of cognitive radio channels 

In between the link from primary to SUs the instantaneous 
channel power gains are represented as g tL (V) , (i0 , 

g 2 i Ci "- 1 and 022 (*0 where v denotes the fading index. By 
using standard signal processing techniques, cooperation of 
PUs, estimation of data or reference signals the CSI of ST-PR, 
PT-PR, PT-SR can be obtained. It is obvious in realistic 
scenarios, that all CSI would contain errors related to 
quantization error, estimation error, out-dated error. We 
assume perfect CSI[8],[9]. 


A. Transmit Power Constraints & Qos Requirement 

Two kinds of constraints such as the peak transmit power limit 
and mean transmit power limit has been used for ST to model 
the power limit for SU’s transmission as 

rJ5(v)< ju*,v v > o (i) 


E{TP 2 (v)}<P« v yv>0 (2) 


Where and P zv represents the peak and mean transmit 
power limits for SU, E ■ denote expectation. Then the 
signal-to-interference plus noise ratio (SINR) outage 
probability constraint for PUs has been adopted 



iv) g 2 1 h :0 






Where and refers the outage probability for PR and 
target SINR, and a p represents the outage probability 

threshold for PR and white Gaussian noise power at PR 
respectively. The PU outage probability when the SU link is 
inactive is expressed as 

f p = Prp 3n M / a? < Yo ] (4) 

Here T refers to the constant transmit power for PT. 


B. Mean Energy Efficiency Maximization 

In this section, the investigation is on the optimal transmission 
strategy for SU to maximize its mean EE. Here, the ratio of 
SU’s mean data rate to its mean power consumption is said to 
be mean EE so as to determine the EE of SU averaged over 
different fading states. The formulation on the SU’s mean EE 
maximization problem as 


PI: h(JP 2 (rj = 




(5) 


Where P c refers the constant circuit power consumption. In 
PI the numerator of h(T P 2 (V) \ is a concave function]! 1 ] and 
its denominator is an affine function[12]. Therefore, 
h(jP 2 CiO ! is quasi-concave for TF Z (V) and PI belongs to the 
quasi-concave programming. 

P2-. = E {Btoa (l 4- - 

q.E{CFP 2 (v) FP C } 


( 6 ) 

Where q is a positive parameter and can be interpreted as a 
pricing factor for SU’s power consumption. As for the 
optimal values of PI and P2, the result is as 
Theorem 1: The optimal mean EE value of PI is obtained 
only if there exists an optimal parameter q* in P2 such that 
T(q*)=0 holds. However, the optimal EE equals g*. 


Proof: \-JjP 2 (t0,! in PI is a concave function over an affine 


function[12]. According to this, we can solve P2 to obtain the 
optimal power allocation policy for PI. To facilitate the 
analysis, an indicator function is introduced for PU’s outage 
event at the fading state v: 




TF-f. uOj? 2 i f lO +.^ 



otherwise 



(3) can be rewritten as 

E{>j, CO }<z s (8) 

Since T(q) is a concave function for, the optimal power can be 
derived since for TP 2 (t0 , T(q) is concave by forming the 
Lagrangian of P2 as 

TP 1 tv)g 12 rv) 

TP l9ll ivU^n 

- q.E[ lTP 2 (v) + P c } -X (E{TP 2 (v)} - -p 

(e {,?p CO} - £ 0 ) (9) 

Where A and // are the positive Lagrangian multipliers with 
(2) and ( 8 ). Then, P2 can be expressed using the Lagrangian 
dual function as 

S = t'lliX T?2 [v]zQ i TP 2 ^ ( 7^2 (v), A >0 

( 10 ) 

Thus, the P2 dual problem is defined as 
P3:min^ 0( ^ 0 g(A f (11) 


L(T P 2 (i-j, A ti) = E | B' log 2 


With a feasible point P2, the condition is satisfied and the 
dual gap is 0 , that is, d* = r* for the convex optimization 
problem. This assures that we can solve P3 to obtain the 
optimal solution to P2. To start with, we consider the problem 
g(k,ju) in (10) for some given X and //. It can be written as 
s (A ii )=E[s '(v) }-)- qP c + /P ; - v + 11 ( 12 ) 


Where, 

9iy) = 

max 7?l Blog^ ( l + 


TP ± B i2 hfi+a£J 


Cgr + ZJTPj Cl?) - fiz p Cp) 


(13) 

Therefore, g(X,ju) is derived by solving the sub problem 
g (rj for each v. For various values of v, maximization 
problems (13) have similar structure. For brevity, we 
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eliminate the fading state index v and the indicator function as 
a function of 7P 2 

max f(TP 2 > = 5^ (l +■ ~ <*£+ &TP 2 

(14) 

is concave of TF 2 . Depending upon KKT conditions, 
the optimal transmit power is derived from (14) as 

* = {[( fl L-+^]n 2 - Tzz Fm * x ) (l5) 

With a turning point y, indicator function is a step function of 




9 2 i V Vp 



(16) 


% (TF Z )=1 if y<0 for any FF z >0. TP 2 * is denoted as the 

optimal solution to the problem. The values of 7P 2 * are 
derived based on the values of x and y in various cases. 

The subgradient method is adopted to derive the optimal 
Lagrange multipliers A.*and //*, that is, we iteratively update X 
and p. 

A (ic + L) = [a* - s k (P. le - E{TP 2 (V)})] (17) 

P {k+1L) = - E {jptO})] (18) 

where k denotes the iteration index and s K refers to a small 
positive step size for the k-th iteration. When the step size is 
constant, the sub gradient method is ensures to converge to 
the optimal value. 


C. Algorithm Steps 

1. We do not need to reorder the MIMO sub channel gain 
realization in a descending order 

2. Take the inverse of the channel gains 

3. Water filling has non uniform step structure due to the 
inverse of the channel gain 

4. Initially take the sum of the Total Power and the Inverse of 
the channel gain .It gives the complete area in the water filling 
and inverse power gain 

5. The initial water level by the formula given below is 
decided by taking the average power allocated (average water 
Level) 

6. The power values of each sub channel are calculated by 
subtracting the inverse channel gain of each channel 

7. In case the Power allocated value becomes negative stop 
the iteration process 


III. SIMULATION RESULTS 



Fig.2 shows BER Vs E- Q fN 2 where BER represents Bit Error 

Rate which defines the number of bits per unit time and 
E-q /A r n defines energy per bit to noise power spectral density 

ratio. It is a normalized signal-to-noise ratio (SNR) measure, 
also known as the "SNR per bit". It is especially useful when 
comparing the BER performance of different 
digital modulation schemes without taking bandwidth into 
account. As the description implies, E h is the signal energy 
associated with each user data bit. 



Fig.3 BER Vs SNR 


Fig. 3 depicts the relationship curve between Bit error rate and 
signal to noise ratio. As the SNR increase, BER constantly 
decreases proving the system to be efficient. 
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Fig.4 EE under different P av using Dinkelbach’s algorithm 

Fig.4 represents mean EE vs. AE under various with 
known F^^ = 0.5 W. The EE starts getting increased with 
P^v . But, when P^ v = 200 mW, it remains stable in 
consideration with SU attaining its maximal mean EE when 
F^, =167 mW. When AE decreases below a particular value, 
there exist the tradeoff between SU’s mean EE and AE. All 
others are consistent with more increase of AE , since SU 
achieves its maximal mean EE for the line with = 200 
mW. 
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Fig.5 shows the mean EE vs. A£T under different Fnm . 
Similarly, the EE also increases with AIT. However, it also 
keeps unaltered with the further increase of AE . Besides, 
increase in F^^, can increase SU’s mean EE given P av <P^ V 
where P^ v is the optimal value. Since the EE has an abating 
margin, the increasing amount of EE decreases when F-naz 
increases further. 



Fig.6 EE under different P av using water filling algorithm 
Fig. 6 shows mean EE vs. AE under various F^- using water 
filling algorithm with known Fimje = 0.5 W. The EE increases 
with F^,. But, it remains stable when F- t , = 200 mW, since 
SUs attain its maximal mean EE when P^ v = 167 mW. When 
AE falls below a particular value, there exist the tradeoff 
between SU’s mean EE and AE. All others are consistent PU 
achieves its maximal mean EE for the line with F^. = 200 
mW. It is observed that the energy is more efficient when 
compared to Fig.4. 



Fig.7 illustrates SU’s EE scheme and Rate-Max scheme with 
Pm&x = 0-5 W. The Rate-Max scheme maximizes SU’s mean 
SE and corresponds to F t when q = 0 in particular. As well, 
the mean EE under Rate-Max scheme first gets increased, but 
with constant = 200 mW and drops quickly with F^, = 
500 mW owing to pursuing SE maximization. 


Fig.8 Comparison between Dinkelbach’s and water filling 

algorithm 

Fig.8 depicts that the energy is much more efficient in water 
filling algorithm than Dinkelbach’s algorithm. The range for 
the intervene between mean EE and As diminishes, since SU 
can be limited by (2) or reaches its paramount EE. Finally, the 
mean EE remains constant when further increases because 
it is either restricted by (2) with small FI or reaches its 
maximal value with a greater FI. All channel gains chase the 
exponential distribution since they are observed to be 
Rayleigh fading and the simulation outcome judge the 
performance. 

IV. CONCLUSION 

In case of successive power allocation the number of 
iterations is more, but in proposed water filling Algorithm the 
number of iterations are less. Initial level of the power 
allocated is close to the ideal value, so the results of proposed 
algorithm are better. In case of equal power allocation, power 
is distributed equally among all the sub channels but 
experiences degradation in performance due to the bad 
channels having low value of SNR. In successive water 
filling, the maximum gain can be achieved, if step size is close 
to Zero (very small) but then the iterations would be infinite 
and time taken to allocate the power would be large. The 
results indicate that the proposed water-filling scheme has 
better capacity. Further, the variation of the outage 
probability is given for the system. 
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